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Abstract. We use the stationary phase method to 
determine the path of optical beams which propagate 
through a dielectric block. In the presence of partial 
internal reflection, we recover the geometrical result ob¬ 
tained by using the Snell law. For total internal reflec¬ 
tion, the stationary phase method overreaches the Snell 
law predicting the Goos-Hanchen shift. 
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I. INTRODUCTION 


Optical beams propagating through dielectrics with dimensions greater than the wavelength of light can be 
described by rays obeying a set of geometrical rules. m In this limit, Snell’s law is used to determine the 
relationship between the incidence and refraction angles and, consequently, the paths of optical beams. gHg 
In this paper, we discuss beam propagation through the dielectric block illustrated in Fig. 1 and, as the first 
step in our analysis (Sec. II), we calculate its geometrical path using Snell’s law. 

The analogy between optics mm and quantum mechanics [TOIITT] has been a matter of discussion in recent 
works. The possibility of linking Maxwell’s equations for photon propagation in the presence of dielectric 

blocks with the quantum-mechanical equations for the propagation of electrons in the presence of potential 
steps [141116) allows one to determine, in a simple and intuitive way, the reflection and transmission coefficients 
at the dielectric block interfaces. [iiHsni As the second step in our analysis (Sec. Ill), we obtain the reflection 
and transmission coefficients for s-polarized waves transmitted through the dielectric system of Fig. I. 

Once we obtain the transmission coefficient at the exit interface, we use the stationary phase method (SPM) 
to calculate the optical path by imposing the cancellation (in the oscillatory electric field integral) of sinusoids 
with rapidly varying phase. [2IH23) The calculation of the position of the maximum of the outgoing beam at 
the exit of the dielectric block, based on the SPM (Sec. IV), represents an alternative way to obtain the optical 
path, one that does not require a geometrical analysis. The SPM analysis only requires that we cancel the 
derivative of the outgoing beam phase. However, the use of the SPM as a mathematical tool to obtain the 
path of optical beams propagating into dielectric blocks is not simply a matter of taste. For total internal 
reflection, [531[^ the SPM also predicts the Goos-Hanchen (GH) shift pSIE5] (see Fig. 2). This shows the 
importance of the SPM not only to recover Snell’s law but also to obtain a typical quantum-mechanical effect. 
The SPM, illustrated in this paper for calculating the path of optical beams, is a mathematical tool easily 
extended to other fields of physics in which wave packets play an important role. 

Finally, after discussing our conclusions, we extend our results to p-polarized waves, suggest how to amplify 
the GH shift by building a band of dielectric blocks, and propose further theoretical investigations. 


II. THE OPTICAL PATH VIA SNELL’S LAW 


Let us consider an incoming gaussian beam, with beam waist size wq and wavenumber fc, which moves along 

the z-direction: mm 




Eoe 
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+ y^ \ 

Wq + ^izjk) ’ 


( 1 ) 


where 

V = xe^+yey + ze^. ( 2 ) 

The plane of incidence is chosen to be the yx-plane. The normals to the left/right and up/dwon sides of the 
dielectric block are respectively oriented along the direction of the unit vectors and , see Fig. 1 (a). The 
unit vector forms an angle 9 with e^, 


Gy = By cos 9 + ez sin 9, (3) 

ej = — sin0-I -02 COS0, (4) 

and 02 ^ an angle 7r/4 with bj, 

By^ = ( By -f 02)/'\/2, (5) 

e^. = {-By + Bi)/V2. (6) 

In the new coordinates systems, the components of the position vector r are then related by 

(y* ^ ^ j_ { 1 ^ \ {y \ 1 \ {y \ 

[ J 1/V^y 1/V -sin6» cos6» J [ z J 
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where R represents the rotation matrix. Observe that a; = x = x*. The new coordinates systems, in this 
Section, will be used to calculate the geometrical path by the Snell law and in the following section to determine 
the reflection and transmission coefficients of the optical beam. 

In order to determine the optical path by the Snell law, we first demonstrate that the outgoing beam is 
parallel to the incoming one, then we calculate the exit point and hnally we obtain the distance d between 
the incoming and outgoing beams, see Fig. 1(a). 

When an optical beam falls onto a boundary between two homogeneous media with different refractive 
indices, it is split into two beams. The refracted (transmitted) beam propagates into the second medium and 
the reflected beam propagates back into the first medium. Snell’s law states that the ratio of the sines of the 
angles of incidence and refraction is equal to the reciprocal of the ratio of the refraction indices of the media 
in which the beams propagate. Thus, for the left air/dielectric boundary we have 

sin6* = nsm.'ip > (8) 

where n is the index of refraction for the dielectric material, 9 and il) are the angles that the incoming and 
refracted beams form with e^. In this case, the second medium is optically denser than the first and consequently 
the refracted angle is a real quantity for all incident angles. The beam propagates into the dielectric block at 
an angle ip with the z-axis and an angle tt/A + ip with the z*-axis, see Fig. I (a). Due to the fact that the 
triangles and CP^.^^^P^^ are similar, 

^-Pleft-Pdown = C'-Pright^up = ^ ” V’ = «• (9) 

Consequently, the optical beam forms an angle ip with the normal to the right side of the dielectric block. By 
using the Snell law, we then find that the outgoing beam forms an angle 9 with e^. This implies that the 
outgoing beam is parallel to the incoming one. 

Let us not determine the point in which the outgoing beam leaves the block, i.e. To do it, the better 

coordinates system choice is represented by the system. Without loss of generality, we can choose 

as origin of the coordinates system the incident point on the left (air/dielectric) interface: 

={ 0 . 0 }. ( 10 ) 

By simple geometrical considerations, we immediately get 

’i} 

and 

C = Pa...+ = + (12) 

To obtain the coordinates of Flight need to calculate the intersection between the following straight lines 


^(^upaPiighi) andr(P^^,,C'), 

In the {gj,^ , Gz,} system, the straight lines are represented by 

(^upPlight) : 2* - = tan (J - V') y* - tan (J + V’) , (13) 

^ (Plight^ = - y*- +cV2^ . (14) 

After simple algebraic manipulations, the previous equation can be simplified to 

^(PppPright) : z^=tan(J^-ip'^ y^-bV2, (15) 

r (Plight C) : z* = - y* + c V2. (16) 
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Finally, can be determined by solving the system m and (HU), 



f ib + c)V2 

[ctan(f - 

i’)-b]V2] 

(17) 

[ 1 + tan (f - V’) ’ 

1 + tan 

(f-V') ]■ 

Once known j observing that 

z^ = 

tan ( J - 9 ] 

1 y* 

(18) 


we can immediately calculate the shift between the incoming and outgoing beams. 


tan ( J - 0) y* (Pight) - (bright) ] / ^l + tan^Q-e) 


b+ (b + c) 


sin 0 


\/V? — sin^ 9 


cos 9 — c sin 9 . 


(19) 


Before concluding this section, we note that, depending on the incidence angle 0, the internal reflections can be 
partial or total. Let us explain this in detail by calculating the critical angle which characterizes the limit angle 
between partial and total reflection. As observed at the begin of this section, at the first air/dielectric interface 
the second medium (n > 1) is optically denser than the first one (n = 1) and we always find a refracted beam 
which moves into the dielectric block forming a real angle ip = arcsin( sin0/n) with the 5-axis. At the second 
interface, we have total internal reflection when 


n sm 


K 


> 1. 


( 20 ) 


It is important to note here that for ip > tt/A the refracted beam cannot reach the second interface. This 
represents an additional constraint to be considered in our discussion. By adding this constraint to Eg. (12011 . 
we obtain the condition for total internal reflection, 


. / 1 \ TT TT 

arcsm — —- < ip < —. 

\n J A A 

In terms of the incidence angle 9, the previous condition becomes 

1 — \/n^ — 1 


arcsm 


In Fig. 2(a), we plot the critical angle 


V2 


9r = arcsin 
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( 21 ) 


( 22 ) 


V2 


(23) 


as a function of the refractive index n. This curve separates the partial and total reflection zones. We conclude 
this section, by observing that for 

1 — \/n^ — 1 


V2 


< - 1 , 


which implies 

n > \Ja + 2s/2, 

we always find total internal reflection, see Fig. 2(a). 


(24) 
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III. MAXWELL EQUATIONS AND TRANSMISSION COEEEICIENT 


In this section, by using the Maxwell equations m we calculate the transmission and reflection coefficient 
at each interface. The phase of the outgoing beam will be then used to calculate by the stationary phase 
method the path of the optical beam. 

The plane wave solution of 




Q a + d - — 

^xx ' '^yy ~ ^zz ^2 


is given by 
providing that 


E{v,t) = 0 
{kxX + kyy-\-kzZ — tot)] , 
k = ^kl+kl + kl = w/c. 


(25) 

(26) 


Obviously a convolution of these plane waves is also a solution of Eg. (1^51) . Many lasers emit beams that 
approximate a gaussian profile 


y,2 p 

E{r,t) = Eq / dkxdky exp [— (^kl + ky) Wq/ 4] exp[i {k ■ r — cot)], (27) 

where wq is the beam waist size. Observe that for <C k (paraxial approximation) the previous 

integral can be calculated analytically and the integration leads to Eq. (Ill- 

Due to the fact that the first air/dielectric discontinuity is along the z-axis, it is convenient to rewrite the 
Maxwell equations by using the coordinate system (x, y, z), 


dxx + dyy + dii - r?(z) ^ 


The plane wave solution is now given by 

exp [ i ( kxX + kyy — ujt)] x 

where 


= 0, withn(z) = 


1 for I < 0 , 
n for z > 0 . 


exp [ f fcj z ] + exp [ — f fcj z ] for z < 0 , 
exp[zg 5 z] forz>0, 


(28) 


(29) 


( fci ) " ^ Tz} = [ky ^ sJv?kP- -kl-kl'^ . 

Observe that the y-component of the wave number does not change because the discontinuity is along the 
z-axis m- By matching the function imposing the continuity of Eieft and Ejeft at the point where the 
refractive index is discontinuous, i.e. z = 0, we find 


n ^ kj - qi 
ki + q-z 


and 


T = 

left 


2 /c? 


(30) 


E,...(r,()=0,with„(«.) = { (31) 


ki + qi 

At the second interface it is convenient to use the coordinate system (x, y^,z^:) and solve the Maxwell equation 

dxx + dy^y, + dz,Z, “ ^ 

The plane wave solution is 

exp [i{kxX + qy,y^, - ujt)] x 

where 


exp [iqz,z^]+ exp [ - z qz, z, ] for z^<al\f2, 

Clown exp [ z kz, z^ ] for z^> 


( 32 ) 


( ^^* ) = R-I f ] ^ and { ky, , } = I qy^ , sjE^-kl- q^^ | . 
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As happens for the first interface, the y^-component of the wave number is not modified because the disconti¬ 
nuity of the second interface is along the z*-axis. By matching the function (15^ at z = we find 


D _ 
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V2, 


(33) 


For the up interface, we can use the result obtained for the down interface. By replacing 

(fcz., gz.) - (fez., gz.) 


a a — b 
and —pz —>■ —■;=^ 

V2 


we find the reflection and transmission coefficients for the up interface, 


R.. = 


gz. - fcz 
gz. + fcz 


exp 


2*gz, 


b — a 


and 


rj. ^ ‘2qz, 

gz, + fcz 


■ exp 


- kz,) 


b — a 

'VT. 


(34) 


In a similar way, the reflection and transmission coefficients for the right interface can be directly obtained 
from the coefficients calculated for the left interface. By replacing 

ki O Qi 

and observing that the discontinuity is now located at z = c, we obtain 

2gs 


Qf — kz 

exp[2zyic] 


right 


gz + ks 

Finally, the outgoing transmission coefficient is 


and 


T = 

right 


gs 


■ exp [i (gi - ki) c] . 


T -T R R T '^kjqi f Qz, - kz, 
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ex' 


p[iqz,bV2 +i{qi - ki)c]. 


(35) 


(36) 


The spatial phases of the optical beam in the different regions are 

in) kx X + ky y + kz z = kx X + ky y + kz z , 

left down) kx x + ky y + q^ z = kx x + qy, yt, + qz, z* , 

down -J> up ) kxX + qy^ y. - qz, z. , 

up right) kxX + qy, y. + gz, = kx x + ky y + q^ z , 

out) kx X + ky y + ki z = kx X + ky y + kz z . 

The outgoing beam, which as expected is parallel to the incoming one, is then given by [TSlfT^ 

2 

exp[- (kl + kl) wl/^] X 


£;out(r,t) = 

47r 


4fcs gj 


dfca;d/ty I 7 \9 I I 7 

(gz + kiY \qz, + kz 


gz. - kz 


exp{ i[q^,bV2 + {qs - ki)c + kr - cot]} . 


(37) 


In the next section, by using the SPM, we calculate the position of the maximum of the outgoing beam and 
consequently the position of the optical beam at the exit of our dielectric system. The calculation based on the 
SPM thus represents an alternative way (with respect to the standard Snell law) to obtain the optical path. 
More important, the SPM calculation also allows one to obtain the GH-shift. 


IV. THE OPTICAL PATH VIA THE STATIONARY PHASE METHOD 


The SPM is a basic principle of asymptotic analysis which applies to oscillatory integrals. The main idea of 
the SPM relies on the cancelation of sinusoids with rapidly varying phase. To illustrate this principle let us 
consider the incoming beam given in Eq. (1271) . In order to maximize the integral we impose that 




u (0,0) 




= 0 . 


( 0 , 0 ) 


The derivatives have to be calculated at the maximum value of the convolution function. For the incoming 
optical beam of Eq. dUl), the convolution function is a gaussian distribution, consequently, the maximum of the 
incoming beam is located at 

X = y = 0 . 

This result, obtained without any integration, is confirmed by Eq. (ED- 
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IV.A PARTIAL INTERNAL REFLECTION 


As discussed in section II, for 9 < 6c, we get partial internal reflection. In this case, the outgoing optical beam 
has, with respect to the incoming beam, an additional phase given by 

(j)= qz,b\/2 + {qi - ki)c. (38) 

By using the SPM, in order to maximize the outgoing integral, we have now to impose the following constraints 


A 

dk„ 


{(j) + k ■ r — ujt) 


( 0 , 0 ) 


A 

dL, 


{(p + k ■ r — tut) 


= 0 . 


(0,0) 


Observing that 
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we immediately find 
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_ ky + gi 

qz A . 

dky. 


dk 
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( 0 , 0 ) 
sind 


\/V? — sin^ 9 


+ 1 cos 9 + c 
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\/ n? — 


cos a — c sm ( 
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6 + (6 + c) 


sin0 


'J'n? — sin^ 9 
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(39) 


(40) 


(41) 


Thus, the SPM, recovering the result (I19L represents an alternative way to obtain the geometrical path of 
optical beams. For partial internal reflection, see Fig. 1 (b), the phase (155)1 is the only phase which contributes 
to the SPM calculation, so this shift in the y-coordinate is the real shift seen in an optical experiment. As we 
shall discuss in the next subsection an additional phase appears for total internal reflection {9 > 9c) and an 
additional shift, which cannot be predicted by Snell’s law, has to be considered. 


IV.B TOTAL INTERNAL REFLECTION 

As anticipated in the previous section, for 9 > 9c an additional phase comes from the double internal reflection 
coefficient 


Indeed, by observing that 


gz. - fcz. 
gz. + fcz. 


(42) 


ki = 


1 - 


sin 9 + \J— sin^ 9 
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and using the constraint (1241) . we have <0. Consequently, the additional phase 


if = — 4: arctan 


\M 


has to be included in our calculation. In order to calculate this new contribution, we start from 


dip 


dk 


= -4 


d 


X,y Q‘z^ 4“ 1^2* P dkx^y 


\kzj 


+ \kz,\ 
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qz,—^ - \kzj 
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which, by using the relation 


2 , I, |2 / 2 IN,2 ^ _ qz, dq^. 


ql+ _ i)fc^ 


dkx^y 1 ^ 2*1 dkx^y 


becomes 


dp 4 dq^. 


dkx^y 1 ^ 2*1 dkx^y 

Finally, the additional shift in the y-axis, also known as GH-shift, is 


S = - 


dp 

dk,, 


( 0 . 0 ) 

4 cos 6 ^ sin Q P \fi 


— sin^ ( 


k — 2 +2 sin 9 \/— sin^ o'j (— sin^ 0) 
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( 44 ) 


(45) 


(46) 


(47) 


The SPM allows one to obtain the y-shift of the outgoing beam both for partial and total internal reflection. 


Ay = 


d 

d + 6 


for 9 < 9c \partial internal reflection], 
for 9 > 9c \total internal reflection]. 


(48) 


In concluding this section, we recall that the GH shift {5) cannot be obtained by using Snell’s law. This 
additional shift is due to the presence of evanescent waves in the air zones close to the down and up interfaces 
and cannot be predicted from geometrical laws [32]. It is similar to the quantum delay time discussed in 
quantum mechanics [55H55] . 


V. CONCLUSIONS AND OUTLOOKS 


In this article, we have shown the value of the SPM as a mathematical tool in determining the path of optical 
beams. Our analysis, which is done in section III for s-polarized waves, can be immediately extended to 
p-polarized waves [Tl fTSlfTOllBS] 


^(P) _ 4 kj qi / gz, - 

{qz + n^ki)"^ \qz, + n^kz^ 


exp [i qz, b V2 + (gi — ks)c]. 


(49) 


For partial internal reflections, the SPM reproduces the geometrical result predicted by Snell’s law, i.e. 
Ay = d. For total internal reflections, the SPM takes into account the GH shift predicting Ay = d-\- 5. This 
additional shift is proportional to the wavelength of the incoming beam, see Eq. (SZl). The order of magnitude 
of the GH shift for a double total internal reflections is thus relatively small and this makes any experimental 
observations difficult. Note that red (A « .633 pm) lasers, whose beam waist is wq = 1 mm, undergo a shift 

— 4 

(5 « 10 wq. Due to the fact that this shift depends on the number of internal reflections, it is clear that, to 
make this experimental measurement possible m, we have to amplify this effect by considering, for example, 
a band of N dielectric blocks. In this case, the final GH shift will be given by NS. 
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To guarantee two internal reflections in each block, we impose that the z,-component of the exit point, 
-^rightJ be the same as the incoming one, This implies, see Eq. uni), 

c tan[^ — tp] = b 


which after simple algebraic manipulations leads to 


c 


T? + 2 sin 0 \/V? — sin^ 0 


b. 




(50) 


In a forthcoming work, we intend to analytically solve the integral (I37p and obtain the outgoing beam 
profile. This analytical study can be done by approximating the transmission coefficient in view of the result 
obtained in section IV. 
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SNELL LAW AND OPTICAL PATH 



(b c) sin 26 



Fig. 1 


Figure 1: Geometric layout of the dielectric block used to determine the optical path by the Snell law. The ? 
and 2 * axes represent, respectively, the normal to the left/right and up/down interfaces. The origin is chosen 
at the point in which the incoming beam touches the first interface, In (a), the 0 axis is obtained from 

the z axis by a clockwise rotation of angle 9. In (b), the rotation is anticlockwise (6* —>■ —9). 
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TOTAL INTERNAL REFLECTION 


CRITICAL ANGLE vs. REFRACTIVE INDEX 


TT 




Figure 2: In (a), the critical angle, dc, is plotted as a function of the refractive index, n. The forbidden region 
represents incidence angles for which the refracted beam at the left interface cannot reach the down boundary. 
For n > v4-|-T\/2, we always find total internal reflection. In (b), we can see that, for 9 > 9c, the geometrical 
path predicted by the Snell law suffers an additional shift, <5. This additional shift, known as Goos-Hanchen 
shift, can be calculated by using the stationary phase method. 
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